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1.0 


INTRODUCTION 


The  possible  utility  of  storage  rings  as  electron  sources 
for  efficient  high  power  free  electron  lasers  has  been  under  study 
for  a  number  of  years.  Studies  of  the  storage  ring  behavior  of 
uniform  one  dimensional  wlgglers  of  the  sort  employed  in  the 
original  Stanford  free  electron  laser  showed  them  to  be  unsuitable 
for  this  purpose.^ The  limitation  arises  from  the  fact  that, 
while  electrons  passing  through  the  wiggler  on  the  average  deliver 
energy  to  the  laser  field.  Individual  electrons  may  gain  or  lose 
energy.  As  a  result,  emerging  electrons  have  a  larger  energy  spread 
than  entering  electrons.  The  growth  of  energy  spread  has  to  be 
damped  by  the  synchrotron  radiation,  a  circumstance  which  limits  the 
laser  power  output  to  be  a  small  fraction  of  the  synchrotron 
radiation  and  hence  limiting  the  overall  efficiency. 

Attempts  were  made  to  remedy  this  situation  by  introducing 
axial  variations  in  the  wiggler  structure,  but  these  proved  to  be 
unsuccessful.  An  understanding  of  this  lack  of  success  was  provided 
by  the  gain-spread  theorem  discovered  by  John  Madey^  and  proved  with 
increasing  generality^"®  over  the  years.  This  theorem  showed 
induced  energy  spread  to  be  an  intrinsic  part  of  any  one  dimensional 
FEL  system  operating  close  to  a  linear  regime  and  hence  subject  to 
the  efficiency  limitation  noted  above. 


The  gain  expanded  free  electron  laser  was  proposed  as  a 
solution  to  the  above  described  problems.  In  its  originally 
proposed  form  it  consisted  of  a  standard  transverse  FEL  wiggler 
modified  so  that  the  wiggler  field  acquired  a  transverse  gradient 
which  Included  a  non-alternating  component.  Electrons  of  different 
energies  were  to  traverse  the  wiggler  in  different  transverse 
positions,  with  the  transverse  gradients  arranged  so  as  to  provide 
equal  transit  velocities  for  all  of  the  electrons.  Since  the  FEL 
amplification  process  is  based  upon  a  velocity  resonance,  it  was 
expected  that  gain  could,  in  this  way,  be  made  energy  Independent. 
Since  the  transverse  gradient  wlgglers  are  two  dimensional 
structures,  the  proofs  of  the  gain-spread  theorem  did  not  apply  to 
them,  and  the  insensitivity  of  gain  to  electron  energy  was  thought 
to  make  efficient  storage  ring  operation  possible.^  Subsequent 
analysis  of  the  uniform  transverse  gradient  wiggler  in  the  form 
originally  proposed  showed  that  while  energy  independent  gain  could 
be  achieved,  the  amplification  process  acted  as  a  driver  for 
transverse  betatron  oscillations.^*^  Numerical  analysis  suggested 
that  this  would  result  in  limitations  for  storage  ring  operation 
similar  to  those  of  one  dimensional  wlgglers.  There  followed  an 
extended  series  of  Investigations  of  various  two  dimensional 
configurations,  which,  while  incorporating  the  principle  of  gain 
expansion,  attempted  to  avoid  transverse  excitation  (i.e.,  the 
excitation  of  betatron  oscillations).  The  systems  were  complex  and 


2 


the  analyses  ultimately  numerical,  often  containing  plausible  but 
not  well  controlled  approximations.  Some  of  them  produced  quite 
encouraging  results,  but  there  was  no  plausible  pattern  of  behavior 
with  respect  to  parameter  variation,  and  no  physical  understanding 
of  what  purported  to  be  successful  designs  was  provided  by  the 
numerical  analyses. 

In  a  recent  report, ^  Rosenbluth  and  Wong  obtained  exact 
analytic  results  for  one  of  the  configurations  that  had  been  studied 
numerically.  These  results  Included  exact  gain  formulas,  Manley 
Rowe  relations,  and  an  extension  of  the  gain-spread  theorem  to 
Include  transverse  excitations.  These  results  proved  to  be  of 
considerable  assistance  In  refining  the  numerical  analysis  that 
preceded  tl  i. 

Previous  unreported  analyses  as  well  as  the  above  stated 
result  suggested  that  the  galn-spread-excltatlon  theorem  for  two 
dimensional  systems  was  of  comparable  generality  to  the  one 
dimensional  gain-spread  theorem.  In  sections  2,  3,  and  U,  we  shall 
show  this  to  be  the  case  and  will  see  that  It  (1)  Implies  the 
Impossibility  of  designing  a  two-dimensional  wlggler  which  produces 
linear  gain  simultaneously  with  vanishing  lowest  order  spread  and 
excitation,  and  (2)  a  wlggler  which  produces  vanishing  lowest  order 
excitation  Is  governed  by  the  gain-spread  theorem.  In  section  5  we 


shall  discuss  the  Implications  of  these  results  for  the  quasl-llnear 
operation  of  storage  ring  FELs  and  conclude  that  they  constitute  a 
possibly  Insuperable  obstacle  to  the  avoidance  of  the  limitations 
found  In  refs.  1,  2,  and  3. 


2.0  THEORETICAL  FRAMEWORK  AND  STATEMENT  OF  THE  GAIN-SPREAD- 

EXCITATION  THEOREM 

As  in  ref.  6  we  begin  with  a  Hamiltonian  of  general  form 

K(E,t,p,x;z)  =  K  (E,t,p,x;z)  +  K  (E,t,p,x;z)  +  .  .  . 

°  '  <2 

where  z  is  taken  to  be  the  independent  variable,  and  Hamilton's 

Equations  take  the  form 


dE  ^  3K  dt  ^  ^ 

Hz’  Tt  Hz’  3E 


^  _  3K  dx  _  3K 

dz  Tx  d?  ^ 


(2.2) 


The  time-independent  Hamiltonian  describes  the  motion 

of  electrons  in  the  wiggler  in  the  absence  of  radiation,  and  z  is 
the  axial  coordinate  of  the  wiggler.  While  we  have  magnetic  fields 
primarily  in  mind,  electrostatic  field  may  also  be  present.  The 
treatment  will  be  sufficiently  general  that  could  include  the 

static  storage  ring  fields  as  well  as  the  wiggler  fields  and  z  can 
represent  any  conveniently  chosen  variable  (e.g.,  the  azimuthal 
angle  about  the  center  of  the  storage  ring)  that  describes  the 
progress  of  the  electrons  around  the  ring. 


terms 


Radiation  is  taken  into  account  by  adding  time-dependent 
Y,  to  the  Hamiltonian  of  order  e”  ,  where  e  is  some  small 
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{parameter.  The  time  averages  of  the  are  assumed  to  vanish,  and 

some  additional  conditions  on  the  time  dependence  will  be  specified 
later . 


We  designate  a  one  parameter  family  of  orbits  (the 
parameter  is  E  )  of  Kq  as  nominal  orbits,  represented  by 
Xq(z,E)  ,  Pq(2,E)  ,  and  tQ(z,E)  with  tQ(0,E)  equal  to  zero.  If 

one  were  considering  the  entire  storage  ring,  we  would  specify  these 
orbits  to  repeat  themselves  on  performing  a  complete  cycle  around 
the  ring,  so  that  they  would  be  periodic  in  z  with  period 
corresponding  to  the  complete  cycle.  Alternatively,  if  the  wiggler 
structure  were  periodic,  one  might  choose  them  to  have  the  period  of 
the  structure.  In  an  ideal  gain  expanded  design  one  arranges  to 
have  tg  Independent  of  E  .  Since,  however,  we  wish  to  be 
general,  we  do  not  specify  any  additional  conditions. 

While  E  is  a  constant  of  the  motion  in  the  absence  of 
radiation,  the  presence  of  radiation  will  cause  it  to  change.  The 
aim  of  a  gain  expanded  wiggler  design  is  to  arrange  things  so  that  a 
particle  entering  on  a  nominal  orbit  of  energy  E  and  changing  its 
energy  to  E'  on  traversal  of  the  wiggler  emerges  on  the  nominal 
orbit  of  energy  E'  .  To  assess  the  extent  to  which  the  objective 
is  achieved,  one  considers  deviations  from  the  nominal  orbit 
designated  by 
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‘  Pg  “  P  ■  Po  •  ^2.3) 

In  the  absence  of  radiation  and  for  small  amplitudes,  a  Courant- 
Snyder  invariant  for  these  variables  can  be  defined  which  is  z- 
independent  and  characterizes  the  extent  to  which  the  orbit  deviates 
from  the  nominal  orbit.  We  refer  to  this  quantity  as  an  action  J 
because  we  shall  construct  a  canonical  transformation  which  reduces 
the  Hamiltonian  to  angle  action  form.  In  the  presence  of  radiation, 
both  J  and  E  became  z-dependent.  Designating  these  z-dependent 
corrections  by  ^  »  where  n  refers  to  the  order  in  e  ,  we 

shall  show  that 

^*^2^  “  7  11“  T'S^ 

o  o 

^*^2^  °  7  1^  7  11“ 

o  o 

We  refer  to  these  two  equations  as  the  gain-spread-excitation  (GSE) 
relations.  The  indicated  averages  are  taken  over  time  and  the 
action  phase. 


3.0  TRANSFORMATION  TO  ANGLE  ACTION  VARIABLES 

Following  the  pattern  of  ref.  11,  we  shall  perform  a  series 
of  canonical  transformations  all  based  upon  a  generating  function  of 
the  form  G(E, t' ,p' ,x;z)  with 


1£ 

3p' 


P 


3G 

ax 


E' 


1£ 

at' 


t 


aE 


K'  =  K  +  ac/az 


The  transformations  to  be  employed  affect  the  K  , 

n 

n  >  1  ,  only  by  virtue  of  the  fact  that  the  variables  upon  which 
they  depend  must  be  expressed  in  terms  of  the  final  canonical 
variables.  We  shall  therefore  keep  track  of  these  relations  among 
variables  but  retain  only  in  the  discussion  given  below. 


Recall  that  Xg(z,E)  ,  Pq(z,E)  , 

conveniently  chosen  one  parameter  set  of  solutions  to  the  equations 


of  motion  generated  by  .  We  begin  with 


(3.2) 


which  yields 


Pg  =  P  '  Po 


(3.3) 


X-  =  X  -  X 

6  o 


(3. A) 


3p„  3x„ 

t  -  t  +  V  °  -  r>  2. 

^  ^  ^  sr-  Pg  3E" 


(3.5) 


E^  =  E 


(3.6) 


3p  3x 

^oa  =  ^o^^'Po  Pg’^o  3l^  ~  Pg 


(3.7) 


In  order  to  simplify  (3.7),  we  make  a  small  amplitude  expansion  of 
Kq  about  the  nominal  orbit.  Hence  we  write 


3K  3K 

.  r  ^  oo  oo 

K  1E,P  +  P„,X  +  X-,z)  =  K  +  X-  -  +  P„  -  + 

o'-  ’*^o  *^6’  o  B  •'  oo  B  7x  ‘^B  3p 


2  2 
,  3  K  3^K 

1  oo  2  oo 

7  -I-  ’‘6  +  — —  Vt  * 


3x 


dx  dp 

O  O 


3^ 

1  oo  2 


*  »  • 


3p  3x 

„  o  o  1  .  2 

=  K  ~  X-  TT -  +  p.  IT -  +  TT-  Ax  + 

oo  B  3z  *^B  3z  2  B 


+  B  XgPg  +  7  C  p2  (3.8) 

where  =  K^(E,p^,x^,z)  .  Although  we  have  omitted  higher  order 

terms  in  x  ,  p  in  (3.8),  we  shall  not  do  so  in  what  follows. 

P  p 

This  is  justified  because  the  higher  order  terms  may  vanish  or  make 
smaller  contributions  than  higher  order  terms  which  will  appear 
subsequently . 1 1  For  compactness  we  have  designated  the  second 
derivatives  by  A  ,  B  ,  C  .  They  are  to  be  thought  of  as 
determined  functions  of  E  and  z  .  We  assume  A  to  be  non¬ 
vanishing  and  positive.  To  complete  the  first  transformation,  we 
combine  (3.7)  and  (3.8)  and  obtain 
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o  TT- 


•^oa  '  I  A  Pe  ^  »  ’'e  Pfi  T  *^00  * 


Next  we  choose 


E  .  3p  . 


n 


oW^  ‘■o 


which  yields 


Pg  “  Pg 


*g"*g 


-  3p  3x 

“  '^g  *^0  *g  “Sir  ■  Pg  -gr 


Using  the  equations  of  motion,  we  see  that 
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i  -i  . 


(3.9) 


(3.10) 


(3.11) 


If  (’‘o  "  ‘'o) 


8x  3p 
o  c 


(  "o  3z“  j  " 


3x  3p  3t 
o  o  o 

^  7z~  7z- 


3p  ,  3K  3p  3K  3x 
o\  oo  o  oo  o 

■5^  I ’^o  75"  "3^“  3jr—  Tir 


3  /. 


/  3K  \ 


,x 

o  o 


'  Ie^  (’^00  *0  75^  ) 


(3.12) 


Hence 


|£--(k  +x  ^ 
3z  \  00  o  az 


(3.13) 


and  completing  the  second  transformation  we  find 


•^ob  "  Y  ^  Pb  ■*■  ®  ‘’e*B  7  ^  *6 


(3.14) 


This  Hamiltonian  can  be  reduced  to  the  harmonic  oscillator 
form,  kg(P^  +  Q^]/2  ,  by  means  of  a  linear  transformation.  Thus  we 


write 


Gc  =  I  P^fj  +  ^  x|  f2  +  Xg  P  fj+  tE 


1  „2 


(3.15) 


which  yields 
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Q  .  Pf^  -f  X3f3 


Pp  “  ’'6^2  '’^3 


:  ^  1  t.2  ^  1  2  ^^2  ^  „  ^^3 

’'p  “  "^  ■*■  7  ^  ^  T  *p  ■55“  ^p*^  3^ 


E  -  E 
T 


(3.16) 

(3.17) 

(3.18) 


We  determine  the  functions  by  requiring  that  (3.16)  and  (3.17) 

transform  a  solution  Q  =  cos  k^z  of  the  harmonic  oscillator 

Hamiltonian  to  a  solution  of  the  equations  of  motion  of  K  ,  .  To 

ob 

this  end  we  make  use  of  the  easily  verified  fact  that  a  K  . 

ob 

solution  may  be  expressed  in  the  form 


X  =  a  /SF  cos^i 
& 


(3.19) 


where 


z 

'I'  =  / 


dz' 

P 


(3.20) 


and  P  is  a  solution  of 


I  P  P*'  -  ^  P*2  +  AC  p2  .  1 


(3.21) 
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with 


P  -  2B  e 


(3.22) 


In  the  above  expression  a  Is  an  arbitrary  constant  to  be 
determined  along  with  f's  ,  and  the  prime  In  (3.21)  and  (3.22) 


Indicates  differentiation  with  respect  to  z  .  Note  that  solutions 
of  (3.21)  can  never  vanish  and  may  be  taken  positive.  Using 


Hamilton's  equations  to  determine  the  related  P  and 
applying  (3.16)  and  (3.17)  we  find  a  “  1  and 


and 


f^  -  -  tan^  (3.23) 

f2  -(|  B*  -  tant)/A8  (3.24) 

fj  =  sec^/ZT^B"  (3.25) 

<|i  =  i|/  -  k^z  .  (3.26) 


Combining  these  with  (3.16)  and  (3.17)  we  find 


Xg  “  /TtF  (Q  cos(^  +  P  sln<}i)  (3.27) 
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r 


Pb 


cos^  -  sin^ 


)+p(i6* 


sin^  +  cos^ 


(3.28) 


It  is,  of  course,  not  necessary  to  know  the  origin  of  the 
f's  to  carry  out  the  transformation.  Using  only  (3.15),  (3.16), 
(3.17),  (3.20)  and  (3.22)  through  (3.26)  we  find 

3G 

K  -  K  .  +  ^ 

OC  ob  dZ 

1  ,2  2, 

-  7  kj  (P  .  Q  ) 

+  (q  co«*  +  Pslii^)^(-||3(3«'  -  ip  +  ACP  ^  -  1  j  . 
Now  applying  (3.21)  we  find 


K 


OC 


1 

T 


kg(p2  +  q2) 


(3.29) 


Since  Kqc  z-lndependent 


J 


(p2  +  q2) 


(3.30) 


is  a  constant  of  the  motion.  Inverting  (3.27)  and  (3.28)  and 
substituting,  we  find 
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1 1  M  ( 


6*^  \  2  * 


1 


(3.31) 


which  we  recognize  as  the  Courant-Snyder  Invariant.  Evidently,  for 

a  given  J  ,  all  associated  values  of  p  ,  x  must  lie  on  an 

p  p 

ellipse  of  area  2xJ  .  The  orientation  and  aspect  ratio,  of  course, 

depend  upon  z  .  It  should  be  noted  that  a  whole  set  of  Courant- 

Snyder  Invariants  exists  because  we  have  only  required  that  S 

satisfy  (3.21).  For  a  study  of  storage  ring  operation  one  would 

normally  choose  6  to  have  the  same  periodicity  properties  as  the 

nominal  orbits.  If  such  periodic  solutions  do  not  exist,  then  the 

nominal  orbits  are  unstable  and  unsuitable  for  storage  ring 

operation.  One  would,  of  course,  require  of  a  design  that  stable 

nominal  orbits  exist.  For  a  periodic  linear  structure  one  would 

have  the  option,  as  for  the  nominal  orbits,  of  choosing  a  B 

function  of  the  same  periodicity.  The  parameter  k  is  also 

P 

completely  arbitrary  and  has  no  effect  on  (3.21)  nor  on  the  relation 

between  J  and  (x  ,p  )  .  Again  there  is  a  natural  choice  when 
P  P 

dealing  with  a  periodic  situation,  and  one  would  choose  k  so  that 

P 

^  has  the  same  periodicity  as  6  .  In  that  case  k  represents 

B 

the  frequency  (wave  number,  really)  of  betatron  oscillations,  and, 
as  noted  In  ref.  10  has  an  Important  bearing  on  the  amplification 
process.  None  of  the  above  has  any  bearing  on  what  we  are  about  to 
prove,  however,  as  the  theorem  holds  for  every  choice  of  B  . 
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For  our  final  cransformatlon  we  convert  (3.30)  to  angle 
(0)  action  (J)  form  In  the  usual  way.  Thus 


G .  -  T  E  +  J 
d 


(3.32) 


which  yields 


E  =  E  ,  T  =  T 


0  =■  sin  ^  Q//7T  or  Q  =  /2J  sin0 


(3.33) 


=\/2J  -  =  /TJ 


COS0 


(3.34) 


K  .  =  k„J 
od  6 


(3.35) 


Our  final  set  of  variables  (E,  t,  J,  0)  are  related  to 
the  original  set  (E,  t,  p,  x)  as  follows: 


X  “  /2jas  sln(^  +  0)  +  X 


(3.36) 


p  =  /737SF  ^cos((li  +  0)  +  y  B*sin(4i  +  0)j  +  p^ 


(3.37) 
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t 


3(ti  1  2  3  B* 


=  T  +  t 

O 


3P, 


+  j  JAB  sln2(({i  +  9)  ^ET  ^  /IJAB  sln((j)  +  0)  ^ 


V2J/A^  fcos((t,  +  0)  +  ^  B*sin((J>  +  0)  J  ^ 


(3.38) 


Our  final  Hamiltonian  Is  simply 


H  =  k  J  +  K  (E,t,p,x;z)  (3.39) 

®  n=l  " 

where  (3.36),  (3.37),  and  (3.38)  are  to  be  used  to  reexpress  t  , 
p  ,  X  In  terms  of  t  ,  J  ,  and  0  .  Note  that  despite  the 
singularity  In  G(,  where  cos^  vanishes,  the  connections  between 
the  old  and  new  variables  and  between  K  and  H  are  non- 
singular.  As  an  overall  algebraic  check  we  have  verified  that  the 
Lagrange  bracket  relations  are  satisfied. 
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4.0  PROOF  OF  THE  GAIN-SPREAD-EXCITATION  THEOREM 

In  this  section  it  will  be  convenient  to  choose  k,  to  be 

zero.  As  a  consequence  the  zero  order  (In  e  )  values  of  the 

canonical  variables  (E  ,  t  ,  J  ,  0  )  are  all  Independent  of  z 

o  o  o  o 

and  will  be  taken  to  be  the  Initial  values.  In  addition  we  shall  be 
Interested  only  In  first  order  corrections  and  time  averaged  second 
order  corrections  to  these  quantities.  Recalling  that  the  time 

average  of  K2  vanishes,  we  easily  see  that  only  K|  can 
contribute  and  hence  we  retain  only  this  term  In  the  sum  over  n  . 
Finally,  we  note  that  as  a  consequence  of  (3.36),  (3.37),  and  (3.38) 

Kj  Is  a  periodic  function  of  0  and  hence  may  be  represented  by  a 
Fourier  series. 

We  could  proceed  by  assuming  only  a  steady  state  stationary 
property  for  the  time  dependence  as  In  ref.  6,  but  in  the  Interest 
of  simplicity  we  specialize  to  a  single  frequency  (taken  positive) 
and  write 


e  e 


IWT 


H^(E,J,z)  e 


lm8 


+  complex  conjugate 


(4.1) 


Then  from  Hamilton's  equations 


,  iwT  V  tt  Imo  . 

-“lee  )mHe  +c*c. 

^  ID 


(A. 2) 


.  lu)T  \  ,,  ImB  . 

=  lew  e  >  H  e  +  c.c, 
^  m 


(4.3) 


de  IwT  r  ime  . 


(4.4) 


dx  iwt  V  1®® 


(4.5) 


We  define 


I  (e  ,J  ,z)  =  /  H  (E  ,J  ,z'j  dz' 
moo  *  m^  o  o  ^ 

o 


(4.6) 


Then  for  the  first-order  quantities  we  find 


J,  =  -ie  e  y  m  I 
1  ^  m 


e  +  c  .c . 


(4.7) 


V  T 

1,  =  lew  e  )  I 
1  ^  m 


e  +  c .c . 


(4.8) 


®1  =  ^  ®  i  3T-® 

o 


(4.9) 


^“^o  ,  i“®0  ^ 

«  4^e  +  c.c. 

o 


(4.10) 


Note  that  the  time  average  of  all  first-order  quantities  vanishes. 
The  quadratic  averages  which  appear  in  the  CSE  relations  may  now  be 


computed.  Denoting  averages  over  r  and  Q  by  <  >,  we  find 

o  o 


<J^>  =  2  I  I 
1  ^  m 


(4.11) 


<Eh  =  2  u)^e^  I  I 
1  ^  m 


(4.12) 


<E,J,>  =  -2  uie^  y  m  1 
11  ^  n 


(4.13) 


Because  general  relations  have  proved  to  be  useful  in  checking  both 
numerical  and  analytic  treatments  of  specific  realizations,  we  note 
that  (4.7)  and  (4.8)  imply 


3Jj  3Ej 

TT"  ~  7W~ 
o  0 


(4.14) 


Proceeding  to  second  order  we  have 


dJ„  lujT  im0  /  a  a  \ 

-ie  e  °  I  m  e  °(j,  +  E,  ^  +  iu,T,  +  1^01  +  c.c. 


TJ-  '  "llfT-"  ■"“'‘I 
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(4.15) 


loar  im0 


=  ieu>  e 


XDIO  / 


^1  ‘“’l  +  ‘"»l)  "mo 

o  o  / 


(4.16) 


Averaging,  we  find 


21 


22 


/  S  dz 
■*  m 


3 

■STT 


(4.21) 


(4.22) 


Combining  these  with  (4.17),  (4.18),  and  (4.11)  through  (4.13)  we 
have  the  GSE  relations 


^•^2^  °  T  ^*^1^  7  ^  ^^1^1^  (4.23) 

o  o 

^^2^°T  TT  ^^1^  T  7J~  ^“^1^1^  • 


Since  the  case  of  small  excitation  has  been  of  particular  interest, 

we  examine  the  J  +0  limit  of  these  relations.  It  is  clear  from 
o 

the  relations  (3.36)“(3.38)  that  we  have 


I 

m 


L  (E  ,J  ,z) 
m  o’  o 


(4.25) 


where  L  =  L  (E  ,0,zl  is  finite  or  zero.  Hence  in  the  limit  of 
mo  m^  o’  ^ 

small  have 


<J?>  ■  ■’ollhol^ 

<eJ>  =  |L  (4.27) 

1  I  00 1 
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<E,Jl>  -  2e2  -  |I.,„|2) 


(4.28) 


which  yields  at  J  =0 
o 


<-'2>  -  ''(I’-lol'  |'•-I0|') 


(4.29) 


<E2> 


2  2  ^l^ool  ^  2  f  ,T 

“  ^  “  Ye—  ^  “  t|L- 


-  IL  l^l 
10  10  ^ 


(4.30) 


It  Is  therefore  clear  that  if  ^*^2^  vanishes,  that  we  have  the 

a 

usual  gain-spread  relation  since  the  -^j-  <EjJj>  term  vanishes. 

2  o 

Furthermore,  if  <Ej  >  also  vanishes,  then,  as  has  been  noted 
previously, 5  gg  joeg  <E2>  ,  and  there  can  be  no  gain. 


Some  comment  on  the  time  dependence  assumed  in  Eq.  (4.1)  is 
appropriate.  We  could  easily  have  replaced  Hm  by 

I  e  .  The  proof  is  essentially  the  same;  expressions  like 


I  ll^l  (4.12)  became  ^  u)2  ll^l^  »  etc.,  (4.23)  and 

m  m,k  '  ' 

(4.24)  are  unchanged  and  (4.29)  and  (4.30)  are  replaced  by 


<J2> 


(4.31) 
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Of  course,  the  time  averages  must  be  carried  out  over  time 
intervals  which  cause  the  frequency  cross  terms  to  vanish,  and  if 
the  frequencies  are  incommensurate  the  GSE  relations  may  contain 
errors  of  order  l/uiT  where  T  is  the  averaging  time. 

For  application  to  the  next  section  we  note  that  if  k  is 

p 

not  set  equal  to  zero,  equations  (A. 11),  (A. 12),  (A, 13),  (A. 23),  and 

(A.2A)  are  unchanged.  The  expression  for  1^^^  [Eq.  (A. 6)]  acquires 

a  formal  dependence  upon  k  ,  viz. 

p 

I  [e  ,J  ,z)  =  dz'  H  (e  ,J  ,z';k„)exp  imk„  z'  (4.6') 
moo  moo  p 

o 

but  from  (3.36),  (3.37),  and  (3.38)  one  sees  that  H  contains 
^  -  kpZ  only  in  the  combination  0  +  ^  from  which  we  conclude 

H  (e  ,J  ,z';k„)  =  H  (e  ,J  ,z'lexp  -  imk„z'  (A. 33) 

0  0  0-'  o  o  '  0 

and  therefore  that  I™  in  fact  also  independent  of  Vc-  • 

P 
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5.0  IMPLICATIONS  FOR  STORAGE  RING  OPERATION 

In  this  section  we  present  a  discussion  of  the  implications 
of  the  GSE  relations  for  quasi-linear  storage  ring  operation  of  the 
FEL. 

In  the  absence  of  the  laser  fields,  the  combined  storage 

ring  and  wlggler  can  be  regarded  as  a  particular  example  of  a 

storage  ring  with  a  somewhat  unusual  system  of  static  fields.  The 

usual  concepts  and  principles  of  storage  ring  design  apply.  Thus 

z  will  be  specified  to  represent  a  cyclic  variable  such  that 

(x,z)  and  (x,  z+Z)  represent  the  same  physical  point.  Hence  the 

Hamiltonian  will  be  periodic  in  z  with  period  Z  .  The  system 

will  be  designed  so  that  stable  nominal  orbits  of  period  Z  exist 

and  are  uniquely  defined  for  each  energy  lying  within  some  specified 

acceptance  range.  The  natural  and  standard  choice  for  3  is  then 

the  (unique)  periodic  solution  of  (3.21),  which,  on  account  of  the 

stability  assumption,  must  exist.  The  natural  and  standard 

specification  of  k  is 
p 


k 


0 


1  dzV6(z') 


(5.1) 


which  causes  (J)  to  be  periodic  as  well.  The  conventional  "betatron 

number"  v  is  defined  by  v  »  k  Z/Zit  .  Because  we  require  the 

p 
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nominal  orbits  to  be  unique,  v  cannot  be  an  integer,  and  to  avoid 
stability  problems  one  designs  the  ring  so  that  v  is  quasi- 
irrational,  that  is  to  say,  not  expressible  as  the  ratio  of  small 
Integers.  As  a  result,  a  given  electron  with  a  specified  J  and  E 
tends  to  uniformly  sweep  through  all  action  phases  0  in  successive 
transits  around  the  ring.  Accordingly,  in  the  absence  of  laser 
radiation  the  electron  distribution  function  will  be  6 
Independent.  We  have  neglected  the  effect  of  synchrotron  radiation 
and  the  RF  cavity  in  the  above  but  it  is  clear  that  the  conclusion 
will  continue  to  hold  when  they  are  taken  into  account. 


We  can  formalize  the  above  by  means  of  the  Llouville 
Equation,  which  for  the  zero'th  order  Hamiltonian  takes  the  form 


3f 

o 

3^ 


•^0  7H~  “  ° 


(5.2) 


The  general  solution  to  (5.2)  is 


fo  =  F(0  -  k^z,T,J,E] 


(5.3) 


where  F  is  an  arbitrary  function  of  its  arguments  except  for  the 
periodicity  conditions  which  follow  from 


f^(0,T,J,E,z)  -  f^(0  +  2x,t,J,E,z) 


(5. A) 


i. 
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f^[0,  T  -  T(E).  J,  E,  z  +  Z)  =  f^[0,  T,  J,  E,  z]  (5.5) 

where  T(E)  Is  the  one  cycle  transit  time  of  the  nominal  orbits  and 
is  given  by 


Z 

T  =  /  t  dz  .  (5.6) 

o 
o 

Equations  (5.4)  and  (5.5)  are  consequences  of  the  fact  that  the 
arguments  on  the  left  and  right  hand  side  represent  the  same 
physical  space  time  point  [see  (3.36),  (3.37),  and  (3.38)]. 

We  note  that  the  space  time  structure  of  the  micropulses  is 
contained  in  the  dependence  of  F  on  t  .  On  account  of  the  RF 
cavity  and  syncrotron  radiation  the  distribution  function  at  fixed 
z  will  settle  down  to  a  periodic  function  of  j  with  period  T^, 
determined  by  the  cavity  frequency.  We  simulate  this  effect  here  by 
neglecting  the  energy  dependence  of  T(E)  ,  assuming  T(E)  =  T 

c 

and  requiring  that  F  be  periodic  in  x  with  period  T  .  Hence 
the  steady  state  F  must  have  both  the  2-n  periodicity  in  0  and 
the  Z  periodicity  in  z  .  Thus  F  will  be  invariant  to  any  shift 
A0  =  2it(r  +  vs)  where  r  and  s  are  integers.  If  v  can  be 
expressed  in  the  form  m/n  where  m  and  n  are  relatively  prime 
integers,  then  a  periodicity  in  0  of  2n/n  is  implied;  while  if 
V  is  irrational,  F  must  be  0  Independent.  We  have,  however. 
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assumed  that  n  Is  not  small,  and  noting  that  the  synchrotron 
oscillation-synchrotron  radiation  processes  will  cause  some 
additional  smearing  not  included  in  the  above  argument,  we  conclude 
that  must  be  effectively  9  independent. 

We  now  wish  to  take  account  of  the  presence  of  the  laser 
field.  We  recall  that  the  CSE  relations  involve  an  average  over  x 
and  9  .  In  order  to  apply  them  to  a  discussion  of  storage  ring 
operation  we  need  to  establish  the  appropriateness  of  such 
averages.  We  first  note  that  the  x  variation  due  to  the 
micropulse  structure  is  slow  compared  to  that  due  to  the  optical 
frequency.  The  x  averages  are  really  averages  over  optical  phase 
and  are  not  significantly  affected  by  the  micropulse  structure.  We 
consider  the  distribution  function  to  consist  of  a  zero'th,  first, 
and  cond  order  part  in  e  ,  and  we  consider  their  effect  on  first 
and  second  order  energy  transfer  and  transverse  excitation  processes 
which  take  place  as  the  electrons  pass  through  the  wiggler.  The 
first  order  part,  f^  ,  contains  exp  iux  factors.  Due,  however, 
to  the  rapid  energy  dependence  in  exp  IwT  ,  and  the  effects  of  J 
and  E  changes  Induced  by  synchrotron  radiation  and  the  cavity  on 
cycle  times  for  Individual  electrons,  there  can  be  no  phase 

correlation  between  the  uix  dependence  of  fj  and  that  of  the 
laser  field  present  when  the  pulse  is  passing  through  the  wiggler. 
Thus  second  order  effects  which  could  in  principle  be  produced  by 
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the  action  of  on  the  entering  fj  contribute  nothing  on  the 

average.  Time  stationary  terms  do  appear  in  £2  »  ''^'1-ch  may  also  be 
9  dependent.  However,  since  ££  is  already  second  order,  its 
effect  on  changes  in  E  ,  J  must  be  of  higher  order  (the  average 
effect  is  fourth  order).  One  might  ask  whether  the  laser  could 
Induce  a  cumulative  buildup  of  9  dependent  £2  of  the  same  order 
as  f^  .  We  consider  this  to  be  precluded  by  the  quasi-irrational 
character  of  v  ,  which  has  the  same  0  smearing  effect  upon  £2 
as  it  has  upon  fg  • 

We  now  proceed  to  a  discussion  of  a  phase  and  6  averaged 
distribution  function  F(J,E,2)  «  (We  ignore  the  micropulse 
structure  here.)  Designating  the  wiggler  entrance  and  exit  by 
z  =  0  and  L  ,  respectively,  noting  that  wiggler  induced  changes 
in  F  are  second  order  and  hence  small,  and  carrying  out  a  standard 
Fokker-Planck  type  development,  we  find  that 

AFj^i  F(J,E,L)  -  F(J,E,0) 

-  -  4  -  w 

2  2  2 

Y  br  mr  *  y  • 

BE  3J 

(5.7) 
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We  are  not  going  to  write  out  a  complete  Fokker-Planck  equation  for 
the  entire  storage  ring  and  are  instead  going  to  aim  for  an  upper 
limit  estimate  of  laser  efficiency  based  upon  entropy 
considerations.  To  this  end  we  compute  the  laser  Induced  entropy 
change.  Let 


S  H  -  /  F  In  F  dJdE 


(5.8) 


then 


AS  =  -  /  AF  In  F  dJdE  -  /  AFdJdE 


(5.9) 


Conservation  of  probability  requires  J  AF  to  vanish.  To 
obtain  this  result  formally  from  (5.7),  we  make  use  of  the  fact  that 
F  must  vanish  when  either  J  or  E  are  infinite  to  obtain 

j  AFj^  dJdE  .  /  |^f(<E2>  -  j  -  I  - 

"  T  ^^1^  H  ■  r  Sj  ‘^'^|e=o  ^ 

+  /  l^F  ^<J2>  -  -J  -JY-  -  J  — yg  j- 


^  </> 

7  -JT 
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1. 


Applying  the  GSE  relations  (4.23)  and  (4.24)  and  the  small  J 
formulas  (4.26)  and  (4.28)  we  have 


/  AF  dJdE 


--4/{ 


2  aF 
<^1>  ?E 


<EjJj> 


3F  \ 

aj  ) 


dJ 


E=0 


(5.10) 


The  last  line  follows  from  the  fact  that  the  storage  ring  aperture 
will  preclude  non-zero  values  for  ^  ^  E  =  0  ,  as  well  as 

from  the  fact  that  all  laser  induced  processes  will  become 
Ineffective  at  E  •  0  .  Combining  (5.7)  (5.9)  and  (5.10),  we  have 


7J 


(<j,>f) 


1  a 

7 


aE 


[<Ej>F) 


2  2  ' 

“  mj  ■  7  In  ^ 

a  J 


(5.11) 


Integrating  by  parts  and  applying  the  arguments  used  in  connection 
with  (5.10)  to  eliminate  boundary  terms,  we  obtain 


1  <e2>  ^  <E  J  >  ^ 


1  3E  3F 

r  yi  TF 


’  13213  \3F 

■  y  TT  ^-^1^  ~  IW  ^1‘^1>)‘37 


1  3F  3F 

r  yr  ^TT 


Another  application  of  ti  '  GSE  relations  provides  us  with  the  useful 


7  I  (w) 


J>  i^f 


(5.12) 


Next  we  apply  (4.11),  (4.12),  and  (4.13)  to  obtain 


LS^  =  G 


(5.13) 


from  which  we  conclude  AS  >  0  as  one  would  expect. 

Li 


Now  consider  the  laser  energy  generation  per  particle, 
AE  ,  given  by 

Lj 
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(5.14) 


AE  =  -  /  <E  >  F  dJdE 


Applying  (4.24)  and  integrating  by  parts  we  find 


Again  applying  (4.12)  and  (4.13)  we  have 


From  (5.13)  we  see  that  in  order  to  have  zero  AS  we  must  have 

JLi 

either  |l  1^  or  oj  ili-  ~  m  ^  zero  for  all  m  .  If  this  were 
I  in|  9E  dJ 

the  case,  however,  (5.16)  tells  us  that  we  also  have  zero  AE  . 

L 

Hence  we  have  rigorously  shown  that  gain  (or  loss)  requires  entropy 
production.  In  order  to  have  a  steady  state  the  total  entropy 
production  in  a  cycle  must  vanish.  Since  the  only  entropy  sink 
available  is  the  synchrotron  radiation,  we  see  that  laser  gain  must 
be  limited  by  the  synchrotron  radiation. 


We  now  seek  a  lower  bound  on  the  entropy  production 
associated  with  a  specified  AE^^  .  To  simplify  progress  in  this 
direction  we  restrict  our  attention  to  fully  gain  expanded  systems, 

O 

that  is,  systems  for  which  the  energy  variation  of  , 

2 

<EjJj>  ,  and  <Jj>  within  the  storage  ring  aperture  can  be 
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with 


E  =  /  EF„  dE 
o  •'  2 


(5.21) 


A  lower  limit  on  AS^  may  be  specified  in  terms  of  these  quantities 

by  choosing  F2  so  as  to  minimize  (5.19)  subject  to  specified 
0 

values  of  6E  and  E  .  A  simple  calculation  shows  the  proper 

o 

form  to  be  Gaussian,  for  which  aE^  =  6E^  .  Hence 


'  >-L 


(5.22) 


AE  6E 


Again  making  use  of  the  gain  expansion  assumption  we  find 


dF, 


1  ] 

7  /  <EjJj>  -gj-  dJ 


(5.23) 


and  applying  (A. 11),  (4.12),  and  (4.13)  we  have 


AS,  >  AF,  i  y  /  1^ 
L  L  ^  m 

m 


(5.24) 


?  r  7  1 

AE,  =  -  e  f  1 1  i  urn  ^  dJ 
L  '  I  m|  dJ 


(5.25) 


Combining  (5.24)  with  (5.25)  we  obtain 
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AS, 


1  I 

m 


/  dF  <uFj\ 

I  ”  ax’  3E“/ 


dJ/F. 


2AEj 

6E 


(5.26) 


and  hence  that 


AS 


L 


2AE, 
^  “5E" 


(5.27) 


It  Is  possible  to  design  a  system  and  determine  an  F  for 
which  the  lower  limit  Is  actually  achieved.  The  system  must  be 

IO 

dominates  (with  m  >  1  order 

to  get  gain).  F2  must  be  chosen  to  be  Gaussian.  The  extra  term 
then  becomes  negligible  If 


dF, 


aiF, 


or 


„  0)  uJ 

"  iSE  ■  sar 


(5.28) 


We  have,  of  course,  no  reason  to  believe  that  the  solution  to  the 
Fokker-Planck  equation  would  actually  take  the  form  which  we  have 
found  for  the  minimizing  F  .  The  general  form  found  for  F  ,  Fj  ^ 
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and  F2  are  all  quite  plausible  But  we  see  no  reason  to  assume  that 

<J>  =  ,  as  implied  by  (5.28),  will  emerge. 

(1) 

Similar  relations  hold  when  there  is  loss  instead  of 
gain.  In  particular  (5.27)  continues  to  be  valid  when  AE  is 

Li 

replaced  by  its  absolute  value.  To  achieve  the  limit  one  must  then 
take  m  negative  and  replace  it  by  its  absolute  value  in  (5.28). 

We  now  calculate  the  entropy  change  due  to  the  synchrotron 
radiation.  For  this  purpose  we  may  ignore  the  presence  of  the  laser 
radiation  and  Introduce  the  radiation  damping  of  synchrotron 
oscillations  in  a  phenomenological  way.  We  study  the  distribution 
function  at  z  =  1  a  point  just  before  the  electrons  enter  the 

microwave  cavity,  and  consider  its  change  as  it  makes  one  turn 
around  the  storage  ring.  That  is,  we  study 

=  F(j,E,n,z^  +  Z)  -  F(j,E,n,zJ  . 

(5.29) 

In  (5.28)  we  have  replaced  the  time  by  the  synchrotron 

oscillation  phase  W  ,  and  have  on  the  basis  of  arguments  given  at 

the  beginning  of  this  section  assumed  that  there  is  no  dependence  on 

the  action  phase  0  of  the  betatron  oscillations.  An  adequate 

phenomenological  Fokker- Planck  expression  for  AF  is 

R 
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(5.30) 


""r  •  -  Ie  “  ”  ■  W  -  TO 


wi  th 


AE  =  -  2a„(e-e  )  -  xa 

S'- 

AJ2  =  Vi[e-E^] 


AJ  =  -2AjJ 


Here  AE  ,  Afl  ,  and  AJ  represent  the  change  in  the  variables 
E  ,  ft  ,  J  which  occur  as  a  result  of  passage  through  the 
microwave  cavity,  wiggler  magnet,  and  transport  around  the  ring; 

X  and  V  are  constant  parameters  which  determine  the  synchrotron 
oscillation  frequency,  is  the  equilibrium  energy,  and  Ag  ,  Aj 
represent  synchrotron  oscillation  and  betatron  oscillation  damping 
constants  per  turn.  Noting  that  J  F^dEdJdil  =  0  we  find 


AS  -  -  J  AF„  In  FdJdEdft 
R  •'  R 


-  -2(Ag  +  Aj)  .  (5.31) 
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Now  applying  the  Robinson  theorem^^  plane  orbits  we  have 


(5.32) 


where  Rg  is  the  synchrotron  radiation  per  turn  per  particle. 
Combining  (5.32)  with  (5.27)  we  obtain  our  principal  result 


o 

This  limit  is  essentially  the  same  as  that  found  in  ref.  1. 
for  the  standard  one  dimensional  uniform  wiggler  and  can  be  shown  to 
hold  quite  generally  for  any  one  dimensional  wiggler  which  satisfies 
the  gain-spread  theorem.  Because  a  necessary  connection  between 
gain  and  entropy  production  has  been  established  by  equations  (5.13) 
and  (5.16)  and  because  (5.33)  has  been  established  in  the  two 
extreme  cases  of  no  transverse  excitation  and  full  gain  expansion, 
it  seems  very  likely  that  a  relation  differing  from  (5.33)  only  by  a 
numerical  factor  of  order  one  holds  quite  generally. 


A  question  which  one  might  raise  in  connection  with  this 

result  and  which  will  be  addressed  in  a  future  publication  is 

whether  it  would  be  practical  to  improve  storage  ring  laser 

efficiency  by  increasing  6E/E  from  the  conventional  0.01  usually 

o 
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assumed  to  say  0.1.  Should  this  prove  to  be  feasible  from  the 
standpoint  of  storage  ring  design,  the  question  of  whether 
sufficient  laser  gain  is  available  under  these  circumstances  to 
permit  laser  operations  would  remain.  Transverse  gradient  wlgglers, 
because  of  their  large  bandwidth  in  energy,  may  have  an  important 
advantage  here.  In  the  conventional  zero  gradient  wlggler  the  gain 
tends  to  vary  as  the  square  of  the  reciprocal  energy  aperture,  while 
for  the  transverse  gradient  wiggler  designed  to  operate  at  the 
m  ■  0  resonance  it  falls  off  only  linearly.  Gain  at  the  m  =  1 
resonance  is  not  directly  affected  by  the  large  aperture,  but  one 
must  take  account  of  the  fact  that  to  achieve  the  limit  (5.33)  the 
steady  state  transverse  excitation  must  be  connected  to  the  energy 
aperture  as  Implied  by  (3.28).  Taking  into  account  the  bandwidth  in 
transverse  excitation  one  again  finds  that  the  gain  tends  to  vary 
Inversely  as  the  energy  aperture.  Investigation  to  determine  which 
of  the  available  options  may  be  most  advantageous  is  in  progress. 


6.0  SUMMARY  AND  CONCLUSIONS 

In  this  paper  we  have  proposed  an  extension  of  the  Madey 

gain-spread  theorem  to  two-dimensional  wlgglers  and  shown  it  to  be 

quite  generally  valid.  It  has  the  Important  consequence  that  an  ,FEL 

wiggler  which  yields  gain  must  at  the  same  time  generate  either 

■  ^ 

energy  spread  or  transverse  excitation.  Furthermore  -wA-have  found 
that  in  an  FEL  operating  quasi  linearly  in  a  storage  ring,  that 
laser  gain  guarantees  the  production  of  entropy  with  every  pass 
through  the  wiggler.  Consequently  the  laser  radiation  generated  is 
restricted  to  be  some  small  fraction  of  the  synchrotron  radiation 
and  our  analysis  suggests  that  this  fraction  is  of  the  order  of  the 

'  f 

fractional  energy  aperture.  This  leads  us  to  conclude  that  the 
achievement  of  high  efficiency  steady  state  storage  ring  operation 
in  a  quasi  linear  regime  is  not  possible.  Nonetheless,  it  appears 
that  gain  expanded  wlgglers  may  have  superior  linear  gain  as 
compared  with  conventional  wlgglers  and  thus  prove  to  be  useful  for 
application  in  storage  rings  with  large  fractional  energy  aperture. 
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